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Abstract. The normal form for a system of ode's is constructed from its polynomial 
symmetries of the linear part of the system, which is assumed to be semi-simple. The 
symmetries are shown to have a simple structure such as invariant function times 
symmetries of degree one called basic symmetries. We also show that the set of 
symmetries naturally forms an infinite dimensional Lie algebra graded by the degree 
of invariant polynomials. This implies that if this algebra is non-commutative then 
the method of multiple scales with more than two scaling variables fails to apply. 
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The method of normal form is one of the most important tools in the study of 
nonlinear differential equations [1]. It classifies the vector field near the critical 
point, and the normal form consists of the resonant vectors characterized by the 
'"q ■ linear part of the vector field. 

Here we study the normal form on where the linear part of the vector field 
is assumed to be semi-simple (diagonalizable) for a simplicity (the general case will 
be discussed elsewhere). The main purpose of this letter is to give a method to 
determine the normal form in terms of the symmetries of the linear part of the 
system. We also show that for certain cases the set of the symmetries forms an 
infinite dimensional Lie algebra graded by the degree of invariant polynomials. As 
a consequence of this algebra structure, we see that the method of multiple scales 
with more than two variables fails in general. 

2. 

Let us consider a system of equations for x E C N , 

(IT 

- = X F (x) = F(x). (1) 
where the vector field Xp is given by 

N 8 rl 

The components of Xp, denoted by the lower-case letters {fi(x)}(L 1 , are assumed 
to be 

fi(x) e C°°(0) (3) 
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for some neighborhood O of a critical point xq. We assume x = 0, i.e., X Q (0) = 0. 
In order to classify the critical point, we write the system in the series near the 
critical point, 

^ = Ax + F^(x) + --- , (4) 

where F^>(x) is a vector valued function of homogeneous polynomial of degree I. 
Here we also assume that the matrix A e Mjv(C) of the linear part is diagonalizable 
with the eigenvalues Ai, A2, • • • , A^. This case is refered to be semi-simple. Now we 
consider the normal form of (4) near x = 0. In his thesis, Poincare showed: 

Theorem. (Poincare [1]) If there is no eigenvalue X s satisfying 

N 

X s = ^2ni\i (5) 

i=i 

for some ni G Z>o and \n\ := En^ > 2 (non-resonant condition), then there exists 
a formal change of variables x = 4>(y) with 4>(y) being polynomial, such that the 
transformed system for y is just the linear part of (4), 

Thus the nonlinear parts in (4) can be eliminated by the non-resonant condition, 
and the solution of (4) near critical point can be described by that of the linear 
equation. The equation for y is called the normal form of (4) . For the case including 
the resonances, we have: 

Theorem. (Poincdre-Dulac [1]) The normal form of (4) is given by 

iL=Ay + GW(y) + .~, (7) 

where G^(y) is a vector valued function of homogeneous polynomial of degree £, 
and is a solution of 



L A G^(y) = [l N (Ay) ■ -fj- - A^\y) = 0, 



(8) 



where In is the N x N identity matrix. 



Outlines of the proofs of these theorems are as follows: 

Under the change of coordinates, x = <j>(y), the vector field X F becomes Xq 
such that 

X^oX F = X G oX 4) . (9) 

Then expanding F, G and <j> as 

F(y)=Ay + F( 2 \y) + ---, 
G(y) = Ay + G^(y) + --- , 

At _ „. 1 a(2) / \ 1 /1 n\ 
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(9) gives, on the n th degree polynomial, 

L A cf>^=F^(y)-G^(y), (11) 
where F^ n \y) is determined successively, e.g., 
F {2) (y) = FW(y), 

F {3) (y) = ^ (3) (y) + <P {2) (y) ■ ^ (2) (y) - (y) ■ ^ (y) , 

and so on. In the case of A = diag(Ai, . . . , Ajv), the operator La in (8) is 

La = In - A- (12) 

Then we first note that the eigenvectors of La are given by a monomial of degree 
|n| = Enj, 

N 



4 U] \y) = ee]\y? for£ = l,...,iV (13) 

i=l 

with the eigenvalue A", LaV^ = ^l^i ' 



,i=i 



where G Z>o and is an eigenvector of A with eigenvalue Xg, Aeg = Xgeg. If 
the Xg satisfies the resonant condition (5), that is A" = 0, then the vector (13) is 
called the resonant vector. Therefore, if there is no resonance (kerL^ = {0}), then 
one chooses G^ n \y) = 0, while for the case of resonances, we choose G^ n \y) £ 
kerL^ fl F^ n \y), which is the resonant vector (13) of degree n. In the next section, 
we relate the resonant vector (13) to a symmetry of the linear part of the system 
(4). 



3. 

Let us first define the symmetry of the system (1). 

A. 

dx 



Definition. A vector field X$ = S(x) ■ is a symmetry of (1), if it commutes 
with the vector field Xp, i.e. 



[Xs, Xp] '■ = Xs o Xp — X F o Xs 

= {S(X).^P--F( X ).^1\.^ = 0. (15) 
1 dx dx dx 



The comdition (15) implies that the vector function S(x) is a solution of the 
linearized equation of (1). The symmetry can be found directly from the general 
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Proposition 1. Let x(t) be the general solution of (1), 

x(t) = U(t;C 1 ,...,C N ) 
with N arbitrary parameters Ci. Then the symmetries are obtained by 

S i (x) = ^-(t;C 1 ,...,C N ). 



(16) 



(17) 



Proof. 



dt 



dx dSj 



dt dx 
d dx 



= F(x) 
d 



— = —-F = 



dSj 

dx 
dU 



da dt oa 



oa 



dF = dF 

dx 1 dx 



□ 



Thus the symmetry is a generator of the transformation group with N parameters 
of the solution. In (16), some of the parameters may be determined by the level 
of the invariant functions of (1). For example, an integrable Hamiltonian system 
with d degree of freedom has d integrals in involution, and the d parameters in 
the general solution are given as the levels of these integrals. With an invariant 
function of (1), we also have: 

Proposition 2. Let be an invariant function, i.e., ^| = 0, and X$ be a 
symmetry of (1). Then <p(£(x))Xs with an arbitrary scalar function <p of £ is also 
a symmetry of (1). 

The proof is straight-forward. This proposition is useful for constructing higher 
order symmetries as we will see below. 

For the case of linear system, F(x) = Ax with A e Mjv(C), the symmetry is 
just a solution of the system. If we assume the eigenvalues of A to be all distinct, 
i.e., Aj 7^ Xj for % ^ j, then we have: 

Proposition 3. The vectors A n x, forn = 0, . . . , N—l, give N linearly independent 
symmetries of the linear system F(x) = Ax in (1) with the above condition for A. 



Proof. Let be the eigenvectors of A with eigenvalue A^, i.e., Aei 



N 



A is diagonalizable, e/s are linearly independent. The equations ^ a iA 

i=i 



\(.e£. Since 
X e t = 0, 



for i 



, N, give 



Ai 
A 2 



N-i n 



Ao^" 1 



1 A n 



A 



N-l 
N 













-ClN- 



= 0. 



(18) 



This equation has only trivial solution a\ = for all i since the determinant of the 
coefficient matrix is the Vandermonde determinant Yl (\ ~ ^j) which is not zero 

by the assumption. □ 
The (linear) symmetries given by polynomial of degree one are called the basic 
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written as a linear combination of the basic symmetries. Namely, for a matrix with 
distinct eigenvalues, we have 

Span{A n : < n < N - 1} = {B £ M N (C) : [A, B] = 0}. (19) 

In order to find nonlinear symmetry, especially of polynomial type, we look for an 
invariant function of the system and use the Proposition 2, as will be shown below. 

4. 

Hereafter we consider the system (4) with the linear part having distinct eigen- 
values. Note that the set of matrix with distinct eigenvalues is generic in Mjv(C). 
In this case, we first note that the resonant condition (5) leads to an invariant 
function of the linear system: 

N 

Propoisition 4. If there exists a resonance, i.e., 3A S = Yl n i^i f or 

some rii £ 1>q 

i=i 

with \n\ > 2, then a function 

N 

t(x) = *7 1 ]l x ?> ( 2 °) 

i=l 

with the coordinates Xi on which A = diag(X\, . . . , Xn), is an invariant function. 

The proof to show d£(x)/dt = is straight-forward. Here we remark that if 
n s 7^ 0, then (20) gives a monomial of degree \n\ — 1. Among all the invariant 
functions of the form (20) for each independent n £ Z> , the invariant with the 
minimum \n\ is called primitive. From the Proposition 2, we have: 

Proposition 5. The vector fields defined by 

XT = SFM • I := r(x)(A>x) 4 , for{ ( _ ™f Z >- _ j (21) 

are symmetries of the linear system. Here d is the total number of primitive invari- 
ant functions, {ii(x)}f =1 , and in the coordinates where A = diag{\\, . . . , \n), the 
invariant function £(x) is given by 



N (j) N (j) 

i=l i=l 

d N 
j=l i=l 



(22) 



N . . N C) d ( ■) 

where A Sj . = n i K with Y n i > 2 (resonant condition), and Mi = Y m j 

i=l i=l J = l 

with nf' = nf' — 5ij . 

Note that in the symmetries (21) the number d is just the total number of 
resonances with independent n £ Z> . In the formula (21), we remark that if 

7^ 0, X™ (x) gives a polynomial symmetry for any m £ Z> and £. However, 

if ni) = 0, X™(x) is not a polynomial, but a certain linear combination of Xp(x) 
with mj = 1 becomes a polynomial. In the latter case, there is no higher degree 
polynomial symmetries with respect to the index mj (i.e., ^j(x)x Sj is polynomial, 
but ^J l3 x Sj for mj > 1 is not). We now show that any polynomial symmetry can 

u„ — u,. „ l: — „„ — ™u;„„j-; — „f /m\. 
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Lemma 6. If X$ = S(x) ■ is a polynomial symmetry, then there exists an 
invariant function £(x) of (20) and 



N-l 

s(x) = ax) c ^ x ( 23 ) 

i=0 

for some constants Ci, for i = 0, . . . , N — 1. 

Proof. We first note that S(x) is a solution of the linear equation, i.e., 

f t S(x) = AS(x) . (24) 

Since A is diagonalizable, we take A = diag(Ai, . . . , Ajv) without loss of generality. 
Then (24) leads to the equations for the component Sj(x) of S(x), 

d , dijc d * \ / a \ d / \ 
— Sjlx) = — ■ —Sj(x) = (Ax) ■ —Sj(x) 
dt n ' dt dx JK ' v ' dx JK ' 

N q 

= XiXi-^-Sj(x) 

i=l % 

= XjSj(x), for j = 1, . . ., N. 
This implies that if Sj(x) is a polynomial, it is a monomial in the form, 



N ... 

A3) 



i=l 



i (x)=l[x^ . (25) 
Substituting (25) into (24), we obtain 

' N 

J2n ( ? ) \-\ j )s j (x) = 0. 



This is just the resonant condition (6). Namely the polynomial symmetry results 
from the resonant condition. Then, from the Proposition 4, there exists an invariant 
function of the form (20), from which we have 



N 

,0) 



s j( x ) = Cj(x)xj = y Xj 1 f[x^ j Xj . (26) 

Since the set of the basic symmetries {A n x}^=o gives a basis of C^, a symmetry 
S(x) = (0, . . . , 0, Sj(x), 0, . . . , 0)* can be expressed as a linear combination of the 
basic symmetries. This completes the proof. □ 

It should be noted that the symmetry given by (25) with zeros for other compo- 
nent, S(x) = (0, . . . , 0, sj(x), 0, . . . , 0)*, is nothing but a resonant vector of (8) for 
A = diag(Ai, . . . , Ajv)- Therefore, we have the main theorem of this paper on the 
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Theorem 7. The normal form (7) can be expressed as 

t = Ay+ £ £ <TS?(y), (27) 

|to|>1 0<£<N-1 

where S™(y) for m G Z> and £ = 0, 1, 2, . . . , N — 1 are given fry f^ij, are 

d 

constants determined uniquely from the original system (4), and \m\ = Yl m %- 

i=i 

Let us illustrate the theorem by taking a nonlinear system of M 2 near an elliptic 
point as a simple example. Namely we consider the system, 

1(3-^ ;)(3***w»*-- <»> 

The primitive invariant polynomial of the linear part of this sytem is £(xi,X2) = 
x\ + x\i and the symmetries (21) are given by 

S^(x) = (xl + xlr(^j, and Sr(x) = (xl + x 2 2 ) m ^y (29) 

Introducing the complex variable z = x\ + ix2, the normal form (27) can be written 
in the following familiar form [1], 

^ = iz+ ^C m \z\ 2m z (30) 

m — 1 

where C m = Cfi 1 - iC^. 

5. 

We now discuss the algebra generated by the symmetries (21). A set of vector 
fields naturally forms an infinite dimensional Lie algebra. Here we show that the 
symmetry algebra gives an infinite dimensional graded Lie subalgebra of the vector 
fields. To do this, we compute the commutator of X™ and X™ by writing these 
in the coordinates with A = diag(Ai, . . . , Ajv), i.e., 

/ N \ N ~ 

*r=(nzr'j£A<z,A on 

Then it is immediate that we have: 

Theorem 8. The set of the symmetries (21) of the linear system forms an infinite 
dimensional Lie algebra with the commutator, 

[XT, xf] = (j2 M ^ x ™ +m ' ~ (E m ^ x T +m ' , (32) 

where M' i =Y l n^m'j. 

i*T n u„,.„ „„,. — „i ™„„u„ — -tl l„j~; — /ort\. 
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Remarks, a) This algebra is graded by the integer m G Z> of the degree of the 
invariant functions £j(x); Z d -gradation. 

b) Since = (Ax) ■ is the original vector field, gives a center of this 

N 

algebra (recall Yl = 0). 

i=i 

c) The set of X™ forms a classical Virasoro algebra (Witt algebra), i.e. 



w. x o'\ = (E M * - E M ) x « +m ' 



E^K-m,))^-', (33) 



where iV~W := £ 
i=i 

Thus nonlinear system having resonances naturally carries an infinite dimen- 
sional Lie algebra of the symmetries, while linear system having no-resonance has 
only a finite dimensional Abelian algebra of the basic symmetries. 

For the example (28), we have an algebra of Kac-Moody (u(l)) and Virasoro 
(Witt) type as the symmetry algebra, 

[X™ 1 X™']=2(m'-m)X™ +m \ 
[X^,xf] =2m'X™ +m \ 

[AT,aT']=0. (34) 
6. 

As a consequence of the previous result, we here show that the method of multiple 
scales with more than two scaling variables fails in general. 
Let t n be the scaling variables defined as 

t n = e n t, for n = 0,1,2,..., (35) 

where e > is a small scaling parameter for the variable x. Then, changing x to 
ex, (4) becomes, 



dx 



(36) 



Now, expanding x in the power series of e, equivalently in the homogeneous poly- 
nomials, i.e. 

x = X W + ex^ + eV 3 ) + • • • , (37) 

where x^ is the homogeneous polynomial of degree £. Substituting (37) into (36), 
we have at order 1, 

^ - AxV = (38) 

dt V ; 

which is the linear lized equation of (1), and in the case of A having distinct eigen- 
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At order e, we have 

Now as a standard procedure of perturbation method, we eliminate secular terms 
in by choosing an appropriate equation for ^f^-. The secular term is given by 
the homogeneous solution of (39), that is, the symmetry of (38). Thus we choose 
^ to be 

^_ = S^(x^) = FW(x^) n kerL A . (40) 

In a similar manner, we have = S < - n+1 ' ) (x^) E kerLA, and the equation for 
x^ in the original time variable becomes 

, (l) OO rj (1) OO 

= Y." -7T- = Axil) + ( 41 ) 

n=0 n n=l 

This is nothing but the normal form (27) with y = x^K It is however important to 
note that the equations of higher orders, = S( n+1 ^ (x^), are not compatible 

in general , as was shown before, i.e. [X S ( n +i),X S ( m +i)] ^ 0, or equivalently 

Q2 X (1) d 2 x (l) 

T r\, r\, , for n^m. (42) 



dt dt dt dt 

This implies that the method of multiple scales, where the scaling vriables t n are 
considered to be independent, fails in general. Note however that the normal form 
(27) makes sense. 

As a final remark, we mention that a similar result can be obtained for a Hamil- 
tonian system where the symmetry algebra is given by a Poisson algebra over res- 
onant polynomials. 
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